Current-current correlations in the three-band model for two-leg CuO ladders 
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We study current-current correlations in the three-band Hubbard model for two-leg CuO ladders 
using the density-matrix renormalization group method. We find that these correlations decrease 
exponentially with distance for low doping but as a power law for higher doping. Their pattern 
is compatible with the circulating current (CC) phase which Varma has proposed to explain the 
pseudo-gaped metallic phase in underdoped high-temperature superconductors. However, for model 
parameters leading to a realistic ground state in the undoped ladder, the current fluctuations decay 
faster than the d-wave-like pairing correlations in the doped state. Thus we conclude that no phase 
with CC order or dominant CC fluctuations occur in the three-band model of two-leg CuO ladders. 
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Since the discovery of the high-temperature supercon- 
ducting copper oxide compounds, the anomalous behav- 
ior of the contiguous "pscudogap" phase has been consid- 
ered a key to understanding the superconductivity mech- 
anism in these materials. However, the nature of the 
pseudogap transition and its order parameter remained 
a puzzle. Early ^iSR experiments^ on YBa2Cu306+ K 
crystals showed evidences for the onset of spontaneous 
static magnetic fields near what was called the pseudo- 
gap crossover temperature T*(x). In addition, differ- 
ent photocurrents for left- and right-circularly polarized 
photons in angle-resolved photoemission spectroscopy^ 
were reported for r^Si^CaC^Os+d"- Now recent polar- 
ized neutron scattering* and Kerr effect measurements 
on YBa2Cu3C>6+:E provide new evidence that there is a 
pseudogap phase associated with a novel magnetic transi- 
tion. The neutron scattering experiments observe a phase 
characterized by a magnetic order which does not break 
translational symmetry and the polar Kerr studies find 
the phase transition at T* (x) which breaks time-reversal 
symmetry. There is at present no agreement regarding a 
theory which encompasses all of these observations. 

Theoretically, ground states, in which circulating cur- 
rents (CC) form spontaneously and thus break the time- 
reversal symmetry, have been found in several models^ 
but in the more realistic t — J model the current-current 
correlations decrease exponentially fast on a two-leg lad- 
der£ However, Varma has argued^ that the minimal 
model for a CC state is a doped three-band Hubbard- 
type model with one Cu d-orbital, one O p^-orbital, and 
one O p a -orbital per unit cells. Using a mean-field ap- 
proach he has found that a CC ground state is possible in 
this model if the Cu-0 hopping integral t p d is of the same 
order of magnitude as the nearest-neighbor Coulomb in- 
teraction V p d between the Cu and the O orbitals and 
larger than the energy difference A p d between these or- 



bitals. This CC state, which breaks time-reversal symme- 
try but not translational symmetry, is consistent with the 
neutron scattering experiments, but further additions^ to 
the model are required to obtain results compatible with 
the orientation of the moments and the Kerr rotation re- 
sults. Moreover, since the interaction between particles 
is strong, the mean-field approach cannot reliably deter- 
mine if a CC phase really exists in the three-band model. 

Several studies of this model have been carried out 
to check Varma's theory using methods for strongly cor- 
related systems. Power-law current-current correlations 
have been observed in CuO chains £ A related "staggered 
flux" phase but no CC phase has been found in the weak- 
and strong-coupling phase diagram of undoped two-leg 
ladders^ Recently, a phase with dominant orbital cur- 
rent fluctuations has been reported in the weak-coupling 
phase diagram of doped two-leg ladders^ However, ex- 
act diagonalizations of small square clusters 1 ^ show no 
evidence for CC patterns in the ground state. Thus, the 
existence of a CC order or dominant CC fluctuations in 
the three-band model is still an open question. 

In this paper we supplement our previous studies of 
the three-band model for two-leg CuO ladder a 13 ' 14 by an 
analysis of the current-current correlation functions in 
doped systems for various parameters A p d and V p d- The 
hole Hamiltonian for this model is given by 
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where the operators d\ a and p\ a create holes with spin 
a in the Cu ti-orbitals and the O p-orbitals, respectively. 
The geometry of the two-leg CuO ladder is illustrated in 
Fig. [TJ where the rung and leg O sites represent p y and 
p x orbitals, respectively. The first and fifth sums are over 
all nearest-neighbor Cu-0 pairs while the second sum is 
over all nearest-neighbor p x -p y pairs on O sites. The in- 
dex i runs over all Cu sites in the third term and over all 
O sites in the fourth and sixth sums. t p d is the hopping 
integral between nearest-neighbor Cu and O sites (solid 
lines in Fig. [TJ and t pp is the hopping integral between 
nearest-neighbor p x -p v pairs on O sites. We have chosen 
the phases of the orbitals such that the signs of the hop- 
ping matrix elements are constant. With the minus sign 
convention of Eq. fTJ one has t p d > and t pp > 0. Ua 
and U p are the on-site Coulomb energies for Cu and O 
sites, respectively. We will work in units where t p a = 1 
and use the typical values t pp = 0.5, Ud = 8, and U p = 3, 
throughout J£ In this model an undoped CuO-ladder cor- 
responds to a density of one hole per Cu site. The hole 
concentration per Cu atom is x = 1 + y with the doping 
rate y = N/(2L), where N is the number of doped holes 
(N > 0, y > 0) or doped electrons (N < 0,y < 0) in a 
ladder with L x 2 Cu atoms. 

In this work correlation functions are calculated nu- 
merically using the density-matrix renormalization group 
(DMRG) method^I We have used up to m = 3400 
density- matrix eigenstates to build the DMRG basis. 
The DMRG ground-state energy is estimated to be accu- 
rate to parts in I0~ 2 t p d or better for ladders with open 
boundary conditions and up to 40 x 2 Cu sites [corre- 
sponding to a total of 282 sites (Cu or O)]. As in Ref. [l2l 
we have calculated correlations 



Cra(«l,«2) = (jr(h)js(i2)) 
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of various local currents j r between nearest-neighbor O 
sites or nearest-neighbor Cu-0 pairs which are illustrated 
in Fig. [TJ We have not found any long-range ordered cur- 
rent patterns for any set of the model parameters that we 
have investigated. Current-current correlations always 
decay faster than 1// as a function of distance I = \i% — i\ \ 
in the present two-leg ladder system. Despite the absence 
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FIG. 1: Schematic lattice structure of a two-leg CuO ladder. 
j r (r — a, b, c, d, e) denote various local current operators con- 
sidered in this work, i is the rung index. 




FIG. 2: Log- log plot of the correlation function |ci(Z)| versus 
I for an electron-doped (top panel) and a hole-doped (bottom 
panel) ladder with A p d = 3 and V p d = 0. Triangles, squares, 
and circles correspond to four, six, and eight doped particles 
= 5, 7.5, 10 %), respectively. Lines are guides for the eyes. 
The dashed lines have slope —1 and —2. Inset: Semilog plot 
of |ci(i)| for a ladder doped with two particles (\y\ = 2.5%). 



of long-range order we can search for patterns in the sign 
of the correlation functions c rs . Close to the ladder ends 
these signs fluctuate as widely as in small square clus- 
ters^ In the middle of long enough ladders, however, 
the sign of a given function c rs does note change with the 
distance I = \i2 — In that case, the relative phases of 
the current-current correlations for various directions are 
compatible with the translationally invariant CC-pattern 
9j proposed by Varma.— We conclude that two-leg CuO 
ladders have CC-like current fluctuations. 

The various correlation functions c rs show a qualita- 
tively similar dependence on the interaction parameters 
and the hole concentration. Therefore, we will discuss 



only ci 



hereafter. As open boundary conditions 



are used, the correlation functions c\(l) = c aa {i\, 1%) have 
been calculated using distances I = \i\ — i%\ taken about 
the midpoint of the ladders (i.e., the integer part of ' 1 ^' 2 
equals L/2). We only show results for I < L/2 = 20 
which have been obtained in ladders with 40 x 2 Cu 
atoms, so that edge effects are small. 

We first investigate the evolution of the current-current 
correlations upon doping. Some results for c\{l) versus 
I are shown in Fig. [2] for various hole concentrations x. 
Although there are substantial differences between hole- 
doped and electron-doped two-leg ladders in the three- 
band model^ we have found that the current-current 
correlations are qualitatively similar in both cases. In 
systems doped with two holes or two electrons these cor- 
relations decay exponentially with distance (see the inset 
of Fig. (2). This behavior can be seen for all interaction 
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FIG. 3: Correlation function |ci(/)| versus I for an electron- 
doped (top panel) and a hole-doped (bottom panel) ladder 
(\y\ = 10%) with A p d = 2. Crosses, triangles, and circles 
correspond to V p d = 0, 1, and 2, respectively. Lines are guides 
for the eyes. The dashed lines have slope —1 and —2. 



FIG. 4: Correlation function |ci(Z)| versus I for an electron- 
doped (top) and a hole-doped (bottom) ladder (\y\ = 10 %) 
with Vpd = 1. Crosses, triangles, squares, and circles corre- 
spond to A p d — 0, 1, 2, and 3, respectively. Lines are guides 
for the eyes. The dashed lines have slope —1 and —2. 



parameters that we have used. In ladders doped with at 
least four electrons or holes (\y\ > 5 %), however, we have 
found that current-current correlation functions exhibit 
an approximate power-law decay l~ v with 1 < v < 2 for 
/ > 3. The overall magnitude of the correlation function 
|ci(Z)| is larger for six or eight doped particles than for 
four doped particles We therefore conclude that current- 
current fluctuations are enhanced upon doping. 

We next turn to the effect of the Coulomb interaction 
Vpd between nearest-neighbor Cu and O sites. Figure [3] 
shows the current-current correlation function |ci(Z)| ver- 
sus I for V p d = 0,1,2. (A recent ab initio calculation^ 
suggests that V p d ~ 1 — 1.5 is appropriate for cuprates.) 
The doping is |y| = 10 % (eight doped electrons or holes) 
which is close to optimal doping in high-temperature su- 
perconducting cuprates. We see that the results are sim- 
ilar for hole and electron doping and do not substantially 
change as a function of V p d- Although the overall ampli- 
tude of |ci(£)| is slightly reduced by increasing V p d, its 
order of magnitude does not change from V p d = to 2. 

The current-current correlation depends more signifi- 
cantly on the energy difference A p d than on the Coulomb 
repulsion V p d- In Fig.[5J we show |ci(Z)| versus I for sev- 
eral values of A p d- While it has been generally accepted^ 
that Apd = 2 — 3, Varma has proposed^ that the CC pat- 
terns are stabilized only when A p d < 0(t p d). We have 
indeed found that the amplitude of the current-current 
correlations decreases with increasing A p d- For A p d = 3, 
ci(Z)| is an order of magnitude smaller than for A p d = 0. 

Our data show that the overall amplitude of current 
fluctuations (for I < 20) increases with doping, decreases 
markedly with increasing A p d < 3 but is little affected 



by Vpd- However, to understand the long-range behavior 
of power-law correlations, it is necessary to investigate 
the variations of their exponent v. We have estimated v 
by fitting our numerical data for the correlation function 
\ci{l > 2)| to a function Al~" ', where both A and v are fit 
parameters. As an illustration Fig. [3 a) shows two such 
fits: The first one corresponds to a rapid decay (y ps 2) 
of the current-current correlations while the second one 
yields one of the smallest exponent, v « 1.2, that we have 
found. As we use data for short distances / < 20 only 
and the correlation functions oscillate widely, the fitted 
values of v are not quantitatively accurate. Nevertheless, 
we think that the variations of the fitted exponent v give 
a qualitative indication of the variations in the long-range 
behavior of the corresponding correlation functions. 

Figure [Hb) shows the fitted values of 2 — v for the 
current-current correlations c\{l) in the parameter space 
(Apd, Vpd) for several hole concentrations. (We use the 
deviation of the exponent from its value in a Fermi sea, 
2 — v, as a measure of the strength of current-current 
correlations.) Two clear trends can be observed both for 
electron and hole dopings: The current fluctuations de- 
crease faster for low doping \y\ = 5 % than for high dop- 
ing \y\ = 10 % and v increases with the nearest- neighbor 
coupling Vpd- The dependence of v on the energy differ- 
ence Apd is irregular but a large value A p d > 3 results in a 
rapid decay of current-current correlation functions. The 
smallest exponent v « 1.2 is found around A p d = 1 — 2 for 
electron doping. For hole doping, however, the smallest 
exponent v « 1.3 is found for A p d = — 1. 

In our previous work-H we studied the pairing correla- 
tions in the three-band model for two-leg CuO ladders. 
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FIG. 5: (a) Correlation function |ci(Z)| for (A p d,V p d,x) — 
(2,0,0.9) (circle) and (3,2,0.95) (square). The dashed lines 
are the fitted functions Al~ v . (b) Fitted values of v in the 
parameter space (A p d, V p d) for various hole concentrations x. 
The radius of the solid circles is proportional to 2 — v. 



We found that electron- and hole-doped systems exhibit 
d-wave-like power-law pairing correlations. Therefore, 
the three-band model at low-doping {\y\ < 2.5%) and 
the t — J model have similar properties: Power-law pair- 
ing correlations and exponentially decaying current fluc- 
tuations^ At high enough doping, however, both pairing 



and current power-law fluctuations seem to coexist in the 
three band model. Comparing fitted exponents for the 
current and pairing correlation functions in hole doped 
ladders (5 < y < 10 %), we find that pairing correla- 
tions always dominate (i.e., decay significantly slower) 
for A p d > 2 while current correlations dominate only in 
a small region of parameter space (A p d < 2, V p d < 1) at 
the highest doping rate investigated (y = 10 %). Thus 
in the hole-doped three-band model on two-leg ladders, 
there is a region of enhanced and apparently dominant 
current fluctuations in good agreement with the interac- 
tion parameters proposed by Varma^ 

For these parameters, however, we showed in our previ- 
ous studyi^ that the undoped ladder has only very small 
charge and spin gaps (which probably vanish in the limit 
of infinitely long ladders). Moreover, local spin moments 
are not formed on the Cu sites as holes are not localized 
on those sites at any doping, and thus there is no ten- 
dency toward (short-range) antiferromagnetic order be- 
tween the Cu sites. Therefore, in the regime of the three- 
band model, where dominating current fluctuations are 
found in two-leg ladders, the undoped system is a param- 
agnetic metal or small-gap insulator. In the regime of the 
three-band model where undoped ladders are "antiferro- 
magnetic" insulators (see Ref. [HI) current fluctuations 
are not enhanced [y w 2) or decay faster than pairing 
correlations. We conclude that no phase with CC order 
or dominating CC fluctuations occurs in the three-band 
model on a two-leg ladders with realistic parameters for 
cuprate compounds. 
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